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Abstract 

y . The importance of including experimental resonances in construct- 

rH I ing effective inter-cluster interactions has been investigated. For this, 

we first address the question of how to obtain the analytical properties 

of the Jost function in regions of physical interest on the complex k- 

plane when the potential is given in a tabular form. We then employ 

C^ '_ the Marchenko inverse scattering method to construct, numerically, 

phase equivalent local potentials supporting the same bound state(s) 
but having different resonance spectra which affect the off-shell char- 
acteristics of the corresponding scattering amplitudes. This implies 
that the inclusion of the experimental resonances in constructing a 
potential would change its shape, strength, and range which in turn 
would modify the bound and scattering wave functions in the interior 
region. This is expected to have important consequences in calcu- 
lations of transition amplitudes in nuclear reaction theories, which 
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strongly depend on the behaviour of the wave functions at short dis- 
tances. The influence of Supersymmetric Transformations on the po- 
sition and movement of resonances has also been investigated. 

PACS numbers: 03.65.Nk, 12.40.Qq, 21.30.-Fy, 34.20.-b 

1 Introduction 

A central problem of the theory of nuclear reactions is to find the interaction 
potential between the colliding nuclei that can describe a wide range of scat- 
tering data and provide the wave functions required. Such a potential is not 
simply the result of mapping data to a convenient functional form. Rather, it 
has an underlying physical basis in that it is related to the nucleon-nucleon 
interaction as well as to the structure and dynamics of the interacting nuclei. 

For practical reasons, it is desirable to express this potential in terms of 
quite simple analytic and local form with parameters adjusted to fit a set of 
scattering data. The phenomenological optical model potentials [I[] are of this 
class. The most commonly used form is that of a Woods-Saxon shape and 
for light nuclei that of a Gaussian. As these potentials fit the scattering cross 
section well, the relevant scattering wave function is asymptotically correct. 
Such potentials, however, do not guarantee that their off-shell characteristics 
are sufficiently good to describe equally well reactions that depend on the 
behaviour of the wave function in the interior region. Such reactions are, 
for example, the electro-disintegration and photo-disintegration processes 
which depend crucially on the wave function at all distances. 

When a nucleus-nucleus interaction is constructed theoretically directly 
from a nucleon-nucleon potential, it is usually nonlocal and quite compli- 
cated. Such potentials can be obtain, for example, by using the Resonating 
Group Method (RGM) in which the inter-cluster interaction is constructed 
via the Pauli antisymmetrization p| . To study the resulting nonlocalities one 
resorts to a construction of an equivalent local potential (ELP) using various 
methods. In this way certain off-shell characteristics of the interaction can 
be revealed. For instance, the use of Fiedeldey's Wronskian method [0, Q to 
construct an energy dependent ELP for the nucleon-a RGM nonlocal inter- 
action revealed a repulsion in the interaction region [^ which suggests the 
possible existence of resonances. 



Another way of constructing local nucleus-nucleus potentials is the use of 
inverse scattering techniques. The potential in this case is directly related to 
the available information of the scattering phase shifts and of bound states. 
In the inverse scattering method at a specific partial wave i (fixed-£ inver- 
sion) [|, ^ H and in the absence of bound states, the constructed potential 
is unique. However, when bound states are present one may construct an 
infinite number of potentials which are spectrum and phase equivalent at 
all energies. This is achieved by choosing arbitrarily the bound states nor- 
malization constants which are not available from experiments. However, 
as we shall show in this paper, the asymptotic normalization constants de- 
termine the distribution of resonances in the fc-plane and vice versa. This 
emphasizes the importance of taking resonances into account in constructing 
effective interactions which are usually ignored. 

One of the main reasons for omitting the resonances is the lack of ex- 
perimental information on their positions and widths especially for broad 
resonances. In the past this was aggravated by the absence of an exact and 
yet simple method to study the analytical properties of the corresponding 
amplitude which could facilitate their incorporation into the potential con- 
struction. Recently, however, a new method for direct calculation of the Jost 



function in the complex fc-plane, has been developed p, [l^, 0, [T^. Within 
this method, the bound, resonant, and scattering states can be found by lo- 
cating the Jost function zeros in the appropriate domain of the fc-plane. It is 
based on a combination of the complex coordinate rotation with the variable- 
constant method used to replace the Schrodinger equation by an equivalent 
system of linear first-order differential equations. Since this method is both 
exact and quite simple, it could be used, together with the fitting of phase 
shifts, to construct a realistic potential which reproduces resonance poles as 
well. 

Complex rotation of the coordinate requires knowledge of the potential 
off the real axis in the r-plane. This poses no problem when the potential is 
given in an analytic form. However, there are cases, such as the aforemen- 
tioned inverse scattering method, in which the resulting potential is given as 
a tabulated set of values versus radial points. In this paper, we show how the 
Jost function method of Refs. [^ 0, |Tl|, |12| can be extended to deal with such 



potentials, and then apply it to study the resonances and their effects on the 
off-shell characteristics of phase equivalent potentials obtained numerically 
by the Marchenko inversion and in Supersymmetric (SUSY) transformations. 



In Sec. II we will describe our formalism by briefly recalling the Marchenko 
inversion method, the relevant SUSY transformations, and the exact method 
of obtaining resonances and Regge trajectories. Our results are presented in 
Sec. Ill and our conclusions are given in Sec. IV. In the Appendix we present 
some technical details related to the inverse method and the Jost function. 



2 Formalism 

2.1 Inverse scattering method 

In the Marchenko inversion scheme [^, 0, ^ the potential Ve{r), for each 
partial wave i, is obtained from the relation 

V,(r) = -2^i^ (1) 

where the function Ki{r, r') satisfies the Marchenko fundamental integral 
equation 

("OO 

Ke{r, r') + T,{r, r') + / K,{r, s)T,{s, r')ds = . (2) 

Jr 

The kernel ^^(r, r') of this equation is related to the ^'-matrix Si{k), and 
thus to experiment, via 

J^,{r,r') = — / hi^\kr) [1 - S,{k)] hf\kr')dk- ^ A^,h^t\h^r)hf\hy) 

(3) 
where h[ (z) is the Riccati-Hankel function, A''^ is the number of bound 
states, and A^e is the asymptotic normalization constant |^, ^ for the nth 
bound state with energy Ef^""' = —h%\/2^ where ibn lies on the positive 
imaginary axis of the fc-plane. The A^i can be expressed in terms of the 
relevant Jost solution f) \k,r) 
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dr . (4) 



The S'-matrix needed in Eq. (^ can be parametrized using the convenient 
rational ansatz, 

Siik) = n 1^ ^ (5) 



i. e. using an infinite number of simple poles and zeros. In practice the num- 
ber of Qi in the above formula is limited to N which is sufficiently large to 
reproduce the input scattering phase shifts and bound states. The conven- 
tions and a method to evaluate these A^ constants are described in Appendix 
A. The rational parametrization of the S'-matrix reduces the Marchenko in- 
verse scattering procedure to an algebraic problem since the kernel of Eq. 
(H) becomes separable; this can easily be seen if one performs the integration 
in Eq. (^ using the residue method. 

In the absence of bound states the above scheme is unique, i.e., once a 
good fit to the data from (0, oo) is achieved, one and only one potential can 
be obtained. In the presence of bound states, however, the potential is not 
unique as it depends on the choice of the asymptotic normalization constants 
Ane which are not provided by experiment. When Ani are chosen according 
to Eq. (HI) or, equivalently, obtained via the Jost function fi{k), 
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nl 



.. fii-k) • 
'dfe{k)/dk^ 



(6) 

k=ib„ 



the resulting potential is unique and has the shorter range 0. Any other 
choice can lead to an equivalent local potential which reproduces the same 
on-shell data but it has different shape and range. 

From Eq. (P) it is clear that values of the normalization constants are 
determined by the explicit form of the Jost function, if known. The Ane can 
be easily obtained using Eq. (j^) where the Jost function is parametrized 
using also a rational representation 

A way to select the alphas and betas, for a finite number of terms used in 
the parametrization (j^) or (|^), is described in the Appendix. 

The Jost function thus constructed approximates well the exact Jost func- 
tion on a segment of the real axis of the /c-plane. From this one should expect 
that it would be a good approximation off the real axis as well. Indeed, if two 
analytical functions coincide on any arc of a continuous curve, they coincide 
everywhere in the region of their analyticity JlS . 



There are at least two problems concerning such a parametrization. Firstly, 
the function (|^ actually does not coincide with the exact Jost function on 



the real axis but can only be a good numerical approximation at certain 
points. How fast the deviations are growing when one moves away from the 
real axis is not known. Secondly, since instead of an infinite number of terms 
in Eq. (^, we have to truncate the product to a finite number of terms, 
not all Oj correspond to zeros of the exact Jost function. In other words, 
the fitting procedure employed to evaluate these parameters, may generates 
a number of unphysical zeros and poles which could be simply an artifact of 
the truncation. 

2.2 Supersymmetric Transformation 

The inter-cluster interaction in nuclear physics is usually characterized by a 
repulsion at short distances. Alternatively, the potential can be made deep 
enough to sustain deep bound states which simulate the role of the so-called 
Pauli Forbidden States (PFS) appearing in the RGM theory which are not 
physical. This means that the Levinson's theorem, 5o(0) — (5o(oo) = tt, is 
fulfilled for this system. These states, however, pose a problem when the 
potential is intended for used in few-body calculations as they generate a 
scattering amplitude with different off-shell characteristics. Thus one resorts 
to the removal of the PFS using a subtraction technique |T^ that results in 



a shallow potential. A rigorous relation between deep and shallow poten- 



tials, however, has been established by Baye [jT5| by using two successive 
SUSY transformations, the first in order to remove the ground state and the 
second to restore phase equivalence. Alternatively the ground state can be 
removed using the Marchenko inverse scattering method and assigning a zero 
value to the corresponding asymptotic normalization constant. The result- 



ing potential is also unique [Tg] coinciding with the one obtained from SUSY 
transformations. In what follows, the main equations of the supersymmetric 
transformation will be briefly recalled (details can be found, for example, in 
Refs. |T3|,|T3). 



The supersymmetric transformation of a given potential V^ that has 
an undesirable ground state at E = E^ ' = —h'^bl/2n, results in a phase- 
equivalent potential V^ which does not sustain this state. Two consecutive 
t r ansf ormat ions , 

Vi'^ =^ y,« ^ y/^) , (8) 

are required for this purpose. The first one eliminates the ground state of 
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the original potential, but also changes the phase shifts. Phase equivalence 
is then restored by a second SUSY transformation. The first step of the 
transformation (E), 



Vrir) = Vr(r) -2-^lnvl/^^/(.6o,r) , (9) 



dr^ 



requires the normalized wave function "^l^^^ibo^r) of the ground state to be 
removed. The second step is implemented via 

VPir) = V}'\r) - 2^1n [¥^X^bo,r)^'\^bo,r)] , (10) 

where the wave function \1/), (ibo^r) is the solution of the Schrodinger equa- 
tion with the potential V^ (r) at the same energy E = E^ and has the 
asymptotic behaviour 



^f)(^6o,r) 



r— >oo 



exp(6orj 



The well-known 1/r^ behaviour of the transformed potential V^ (r) at short 
distances can easily be deduced from Eq. (|TD|). Indeed, both wave functions 
\1/q^ (ifeo?''") and W^ {ibQ,r) near the origin have the behaviour ~ r"^ with 
some integer m, which makes the second term of this equation proportional 
to ~ l/r^. 

Finally, we note that the SUSY potentials are intrinsically related to in- 
verse scattering theory |]l^ which can be used to generate energy-independent 



shallow potentials from energy-dependent deep potentials |]T8|. For this, 
one has to calculate first the phase shifts at all energies and then employ 
the Marchenko scheme to construct an ^-dependent potential from which a 
unique shallow potential can be obtained either by super transforming twice 
or, as mentioned earlier, by assuming that the corresponding asymptotic 
normalization constant is zero. 

2.3 Exact method for locating resonances 

The method we employ here for locating potential resonances belongs to the 
class of so-called complex energy methods which are based on a rigorous 
definition of resonances, namely, as zeros of the Jost function. Unlike most 
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of the other methods of this class, which involve an expansion of the resonant 
wave function in terms of square-integrable functions, our method is based 
on a direct calculation of the Jost function at complex k by integrating exact 
differential equations equivalent to the Schrodinger equation. 

By locating a complex zero kj. of the Jost function in the fourth quadrant 
of the momentum plane, we can obtain, at the same time, the resonant energy 
Er and width F from the simple relation 

^' - f = t/' ■ 

It is clear that the accuracy of E^. and F is related to the accuracy in calcu- 
lating the Jost function itself. The latter is obtained, in our method, from 
the asymptotic value of a function F^~'{k,r) , 






Mk) =lim Fr\k,r), (11) 

|r|— >oo 



which, together with its partner F^ {k,r), obeys the system of first order 
differential equations 

drFt\k,r) = b_^v,{r)[hf\kr)F^-^\k,r) + h\-\kr)F^-\k,r)\ , 

(12) 

drFt\k,r) = -h_it!lv,{r)[h^^\kr)Ft\k,r) + h[-\kr)Ft\k,r) 

the boundary conditions being 

Ft\k,0) = l. (13) 



The origin of the relation (|Tl]) becomes clear if one notices that the sum of 
the products of the auxiliary functions F^ (k,r) with the Riccati-Hankel 
functions, 

(Pe{k,r) = - [h['^\kr)Fi^\k,r) + hi'\kr)Fi'\k,r)] , (14) 

obeys the Schrodinger equation. The function (pelk, r) is the so-called regular 
solution which vanishes near r = exactly like the Riccati-Bessel function, 

i.e. 

lim(l)i,{k,r)/je{kr) = 1 . (15) 

r— >0 



The Jost function, Eq. (|r^), can be obtained by comparing the asymptotic 
behaviour of the regular solution, 

Mk,r) = ;^l[h^;-\kr)f;{k*) + h[-\kr)Uk)] , (16) 

with Eq. ([I^ expressed in terms of the auxiliary functions F^ (k,r). 

In Ref. [jlO| it has been shown that the limit ([TT| ) exists for all complex k 



for which 

Im A;r > . (17) 

If r is real, the condition (^) is only satisfied for bound and scattering states 
but not for resonances. Therefore, to calculate fe{k) in the fourth quadrant 
we make the complex rotation of the coordinate in Eqs. (|T^), 

r = xexp{ie) , x>0, O<0<^, (18) 

with a sufficiently large 6. 

The above scheme works extremely well in locating bound, scattering, 
and resonant states as well in finding Regge poles and trajectories when the 
potential Ve{r) is central, short range and it is given in analytic form. A 
generalization to non-central, multi-channel, and Coulomb-tailed as well as 
to singular potentials can be found in Refs. [|, |10|, |ll], [1^. However, in 



this work we are concerned with potentials given in a tabular form and the 
question of how to handle such potentials will be discussed when presenting 
our results. 



3 Results 

In the present work the potentials were generated either by inversion or by 
SUSY transformations and hence in both cases they are available in a tabular 
form. To make an analytic continuation of them into the first quadrant of the 
complex r-plane, needed for the complex rotation, we fitted the potentials 
on the real axis by simple analytical forms with adjustable parameters and 
then considered r in these forms as a complex variable. Such an approach to 
analytic continuation is based on a theorem of the complex analysis which 
states that if a function is analytic in a region and vanishes along any arc of a 
continuous curve in this region, then it must vanish identically in this region 
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1^ . The obvious corollary of this theorem is that if two functions coincide 
on a curve, they coincide everywhere in the region of analyticity. Therefore, 
the analytical form which coincides with the tabulated function on the real 
axis should reproduce this function off the real axis as well. The question 
then arises what if the potentials coincide within numerics. In other words, 
we want to know whether small numerical deviations in the potential on the 
real axis generate perceptible deviations of the position of the resonances. 
We have investigated this situation first by assuming the following analytic 
potential 

V{r) = 5 exp [-0.25(r - 3.5)^1 - 8 exp(-0.2r2) (19) 

where the strength parameters are given in MeV and r in fm. The reduced 
mass /i is such that h'^/2ji = 1/2 MeV fm^. The resonances and Regge tra- 



jectories of this potential were investigated in Ref. [T^. We then fitted the 
N points V{ri), i = 1,2, ■ ■ ■ ,N using the ansatz 

V&t{r) = Y^ an exp -6„(r - c„)^ + ^ rf„exp(-enr^) (20) 

n=l n=l 

where the parameters were obtained using the MERLIN minimization pro- 
gram ll^. The minimization was stopped when the least square error on 
120 points and A^i = 5 and N2 = 4: was of the order of 10~^ - a typical 
accuracy in these cases. This means that the fit to the analytic potential was 
between the third and fourth decimal in the whole region. With such a fit, 
all resonances found in Ref. |1^ were recovered within three to five decimal 
points. For comparison, the energies and widths for few of them (the lowest 
resonances in each partial wave) are given in Table |l]. Obviously the accuracy 
can be improved as the fit to the potential improves. One further comment 
is necessary: The form factors used for the fit restrict the use of a rotation 
to only a certain region. In this respect the use of splines is not suitable at 
all as they diverge for all angles of rotation. 

These test calculations show that the Jost function method based on 
the complex rotation of the coordinate, is applicable and retains its effec- 
tiveness even when the potential is given numerically on the real axis. We 
can therefore use it to study the importance of incorporating physical two- 
cluster resonances in constructing a potential. This can be easily investigated 
by constructing ELP's for a specific partial wave via inversion as described 
above and studying the implications of the implanted resonances. For this we 
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use the nucleon-o; potential of Dubovichenko and Dzhazairov-Kakhramanov 



2^ for the i = partial wave 

V{r) = -Voexp{-ar^) (21) 

where Vq = 55.774 MeV and a = 0.292 fm^^. This is a deep potential that 
sustains an unphysical PFS state at —9.73058 MeV. This means that the 
Levinson's theorem, So{0) — 5o(oo) = vr, is fulfilled for this system. At large 
distances the radial wave function decays exponentially, 

uo{r) — ^ Asexp{-bor) , 

-1 /Q 

and the asymptotic normalization constant was found to be Ag = 6.1603 fm ' . 

By varying the asymptotic normalization constant we obtained a set of 
potentials which were fully phase shift and bound state equivalent but have 
different number of potential resonances. These potentials are shown in 
Fig. ^ It is seen that for values of Ag less than the choice given by (0) 
[As = 6.1603), a hump appears in the interaction region which is higher as 
As becomes smaller while at the same time the well becomes deeper. For 
values larger than 6.1603 the potential is also of long range but without a 
hump. In the extreme case oi Ag = the potential becomes repulsive at all 
distances. This means that as Ag ^ resonances are generated and their 
appearance and position depend on the specific choice of Ag. 



Using the ansatz (^OD we fitted these potentials, with A^^i = 5 and N2 = 3, 



via the MERLIN code |T^ the accuracy being again within a fourth decimal 



at all points meaning that corresponding accuracy in reproducing the phase 
shifts was better than 0.0001 of a degree. We employed the analytical rep- 
resentations of these phase-equivalent potentials to locate the zeros of the 
Jost function in the fourth quadrant of the /c-plane. The original potential 
(I2TD, which is also a member of our set of the phase-equivalent potentials, 
does not generate any physical resonances. All the zeros of the Jost func- 
tion, which we found for this potential, are situated below the diagonal of the 
fourth quadrant of the fc-plane and, therefore, represent sub-threshold res- 
onances which are unphysical. The growth of the potential barrier when Ag 
decreases, indicates that some physical resonances should appear. In other 
words, when Ag becomes smaller some of the Jost function zeros should move 
up to the area above the threshold line. When, however, Ag is too small, the 
barrier transforms into a strong repulsive core, and the resonances should 
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disappear. This can be seen in Fig. |^ where we present the distributions of 
the Jost function zeros for three phase-equivalent potentials corresponding 
to As equal to 3, 2 and 1. It is clear that the choice oi Ag determines positions 
of the resonances and vice versa. 

At a first sight one can argue that the zeros of the Jost function practically 
have no effect on the scattering processes because they are far away from the 
real axis. These potentials, however, generate bound and scattering wave 
functions which have a different behaviour in the interior region. The bound 
state wave functions are shown in Fig. § while the scattering wave functions, 
for center of mass energy E = 5 MeV, are plotted in Fig. ^ The nodeless 
wave function for the SUSY potential V^'^\r) is also shown in the latter. Since 
the interior region (within few fm) is of importance in describing various 
nuclear reactions, the existence and distribution of resonances cannot be 
ignored when the reaction observables are calculated. These differences are 
also a source of off shell differences in the corresponding scattering matrices 
which are manifested in three- and four-cluster calculations. 

As another example, we consider the a-a local potential of Buck et al. 
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4e2 
Vir) = -Voexp(-ar^) + — erf(/3r) , (22) 

r 

with Vq = 122.6225 MeV, a = 0.22 fm'^, and f3 = 0.75 fm^^ This po- 
tential sustains, in the i = partial wave, two unphysical deep bound 
states at -72.78 MeV and -22.28 MeV and a resonance at 179.22 keV with 
r = 0.94905 keV. The position of the resonance is quite different from the 
values of 92.12±0.05keV and F = 5.80 eV given by Buck et al. ^). A pos- 
sible reason for such a discrepancy is that in locating these spectral points we 
included the Coulomb tail of the second term of Eq. (B^) in an exact way as 



was proposed in Refs. 0, |T0[ . We are interested about the movement of this 
resonance when we eliminate one of the unphysical bound states using the 
SUSY transformation. To this end, the numerically obtained SUSY potential 
V^"^^ was fitted using the ansatz 

5 5 

Vfit{r) = Y^ anexp(-6„r^)/r^ + ^ c„exp(-c„r^) (23) 

n=l n=l 

with the accuracy of the fit of 450 points being better than 10~^. Using 
this analytic potential, we found a zero of the S-wave Jost function a.t k = 
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0.10383-z0.77289xl0-^fm-i (^.e. at E = 0.11261 - -0.33532 x 10"^ MeV) 

which is significantly different from the position of the S'-wave resonance 
given in Table 0. In view of the above discussion this is not surprising since 
the SUSY transformation drastically changes the asymptotic normalization 
constant for the ground state, to zero. 



4 Conclusion 

We demonstrated that for a potential given numerically the analytic prop- 
erties of the corresponding Jost function in the complex fc-plane can be 
obtained via fitting the potential by any analytic form that allows a complex 
rotation into the first quadrant of the complex r-plane. The scattering ob- 
servables, the bound states, and the potential resonances can be calculated 
with a sufficient accuracy which is improved with an improved fit to the 
potential. 

The phase shifts and therefore the on-shell S-matrix which are extracted 
from experimental data on the real /c-axis, contain information about reso- 
nances in an indirect way. The phase-shifts "feel" the existence of resonances 
only when they are close to the real axis (narrow resonances). The broad 
resonances, however, remain "unnoticed" by the phase-shifts and therefore 
a potential which is based on them, generates an S'-matrix without the cor- 
responding poles. However, even extremely broad resonances affect the be- 
haviour of the physical wave function at short distances. This implies that an 
information on the distribution of resonances can be a clue for making a cor- 
rect choice among very different potentials generating the same phase-shifts 
and the same on-shell S'-matrix. 

Such information can, in principle, be obtained from various inelastic pro- 
cesses. For example, photo-excitation of a nucleus and its subsequent decay 
in two fragments A and B can reveal y4i?-resonances which are not "visible" 
in elastic 745-scattering. When an effective potential is constructed us- 
ing not experimental but theoretical phase-shifts, additional effort to locate 
broad resonances would exclude ambiguities (non-uniqueness) of the effective 
potential. This is the case, for example, with the RGM theory which can 
produce very complicated (nonlocal) potentials. To apply such a potential 
in realistic calculations, one usually calculates RGM phase-shifts and, using 
them, constructs a simple effective (phase-equivalent) potential. 
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In summary, when constructing a realistic effective interaction, it is im- 
portant to tak;e into account physical resonances. In the case of an energy- 
and ^-independent potential, the resonances at all partial wave must be in- 
corporated. This will guarantee that a reliable local potential is obtained 
that generates a transition matrix which has correct off-shell behaviour. 
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A Distribution of the 5— function poles in the 
complex /c— plane 

The ^'-function is defined in terms of the Jost functions by 

S,.) ^ ff (24, 

The fe{k) can be parametrized using the rational form 

/K^) = n^ (25) 

implying that it has N simple poles a.t k = jSi and N simple zeros at k = ai. 
The f{k) is a well defined function in the upper-half fc-plane. The poles /3j 
are therefore situated in the lower fc-plane. In contrast, the a^ can be in 
both planes. Those in the upper half fc-plane are denoted by Tj and those in 
the lower by a^. For the S'^-function we have 

or 

MA:) = np7^np7^n|-f^ n.+n^=n (27) 

Setting ^i = —Pi (Im ^i > 0) we may rewrite Eq. (^) as 

s<w = nf^ni^nf^ (28) 

or adopting the same symbols for all poles and roots 

Si{k) = n l^ , M = 2N (29) 

m=l '^ '^"i 

It is clear that the number of poles A^^, N^, and A^^- for the ^, a, and r poles 
are related by N^ = N^- + N^-. 
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The form ( p9D is the one used to parametrized the S'-function or equiv- 
alently the phase-shifts 6. This can be easily achieved by rewriting (|2^) in 
the form 

or 

M M 

Se{k) + S,{k) J2 (-1)'"^™A;" = 1 + E ^-fc'" (31) 

m=l m=l 

from which the Am are evaluated by choosing M different values Si{ki). The 
poles Oj are then determined as the zeros of the polynomial in the denomi- 
nator of Eq. (ig). 

Once the aj are found, one has to extract the a's and /5's and their 
distribution on the fc-plane. In the absence of bound states we have 

M =N, + N^ 

from which 

M 
N, = N^ = - (32) 

In practice one varies M and the position ki until a good fitting of Se{k) 
at all energies is achieved that fulfills the condition (|32D. In the presence of 
bound states the sorting is quite tricky. For example, in the presence of one 
bound state one has 

M =N^ + N^ 

N^ = 1 

Thus 

^,, = ^, iV, = ^. iV, = l. (33) 

The number of poles in the upper fc-plane A^^^ is given by 
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while the number of poles in the lower A;-plane A'^; by 

Therefore 

The same argumentation can be used for any number of bound states Nr to 
obtain 

N^ = Ni + 2K ■ 

This condition must be satisfied for a correct construction of a potential. 



Table 1: Comparison of the parameters of the first resonances generated 
in each partial wave by the exact potential (19) and by the approximate 
analytical form (20). 





exact potential 


fitted potential 


£ 


Er (MeV) 


r (MeV) 


Er (MeV) 


r (MeV) 





2.252380731 


0.000118256 


2.2520 


0.0001179 


1 


0.807634844 


0.000000110 


0.8082 


0.000000110 


2 


2.384151637 


0.000082862 


2.3843 


0.000082812 


3 


1.009031953 


0.000000046 


1.0095 


0.000000046 


4 


2.676524768 


0.000027824 


2.6768 


0.00002783 



Table 2: Spectral points generated by the a — a potential proposed by Buck 

et al. [23]. 





momentum (fm ^) 


energy 


(MeV) 


£ 


Re (A;) 


Im(A;) 


Re(E) 


r 








2.63671 


-72.62600 











1.56602 


-25.61905 





2 





1.45123 


-22.00095 








0.13098 


-0.00017340 


0.17922 


0.00094905 


2 


0.53895 


-0.061203 


2.99519 


1.37832 


4 


1.07513 


-0.084009 


12.00121 


3.77407 


6 


1.85894 


-0.48994 


33.59171 


38.05717 
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Figure 1: Bound state and phase equivalant potentials for the nucleon-a 
interaction. These potentials generate different resonance spectra. 




Im k (fm ) 

Figure 2: Distribution of the Jost function zeros (filled circles) in the com- 
plex /c-plane for the S'-wave N-a potential for three different values of the 
asymptotic normalization constant A^. The diagonal of the fourth quadrant 
represents the threshold boundary Re -E = 0. 
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Figure 3: Bound state wave functions generated by the potentials of Fig. 1 




Figure 4: Scattering wave functions generated by the potentials of Fig. 1 for 
the center of mass energy E = 5 MeV. The Ag = corresponds to the shallow 
V^'^^ potential that generates a nodeless wave function at short distances. 
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